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ABSTRACT 


A load relief autopilot design to minimize the maximum thrust vector deflection angle 
required to maintain vehicle stability during flight through design winds is analyzed. 

The feedback variables are angle of attack and flight path angle as well as attitude error 
and rate. A rigid body configuration is assumed, and sensor dynamics are not considered. 
Deflection angle and attitude error responses are presented for design winds. It is shown 
that the load relief autopilot analyzed results in a 50 percent decrease in deflection 
requirements for these winds as compared to a conventional autopilot designed to trim 
the vehicle. 
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SUMMARY 

A load relief autopilot designed to minimize the maximum thrust vector deflection 
angle required to maintain vehicle stability during flight through design winds is analyzed. 
The feedback variables are angle of attack and flight path angle as well as attitude error 
and rate. This autopilot is investigated for a vehicle consisting of a 260-inch solid motor 
first stage, SIVB second stage, and extended Voyager payload. A rigid body configuration 
is assumed, and sensor dynamics are not considered. 

Deflection angle and attitude error responses are presented for expected flight winds. 
It is shown that the load relief autopilot analyzed results in a 50 percent decrease in de- 
flection requirements for these flight winds as compared to a conventional autopilot de- 
signed to trim the vehicle. It is also demonstrated that the results should apply to most 
other vehicles, as well as the one simulated. 


INTRODUCTION 

Due to the limited gimbaling capability of solid fuel rocket vehicles, it is extremely 
desirable to minimize the maximum thrust vector deflection angle required to maintain 
vehicle stability while flying through design winds. Stability is defined here as the ability 
to follow a desired flight path or return to this path after a deviation due to some distur- 
bance. Deflection angle requirements can be reduced by the addition of aerodynamic sur- 
faces (fins or canards; see ref. 1) or by autopilot design. Conventional autopilots are 
designed to maintain trimmed flight conditions (i.e. , to maintain the nominal flight 
path during the disturbance). Usually, the deflection angle is not permitted to reach its 
maximum value. 

Borsody and Teren (ref. 2) define an ideal autopilot which allows the deflection angle 
to increase to its maximum value and remain at this value for some finite time. They 


further determined the vehicle stability limit; that is, the minimum value of the maximum 
deflection angle for which the vehicle is stable as defined above. It was demonstrated 
that with a feasible open loop control, the deflection angle requirement for design winds 
could be reduced to 56 percent of the trim requirement. 

A type of load relief autopilot designed to exploit this advantage is investigated in this 
report. The feedback variables include angle of attack and flight path angle in addition to 
the conventional feedback variables of attitude and attitude rate. The approach taken is 
that a maximum allowable magnitude for the deflection angle o max is assumed. The 
minimum value of required for stability is established by reducing until the 
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vehicle becomes unstable for some expected flight wind. The vehicle uses the commanded 
deflection angle whenever this value is less than 6 max> If the commanded deflection 
angle exceeds the limit, 6 max is used. If the system is stable, the commanded deflec- 
tion angle will return to a magnitude less than 6 max after a short time. 

The analysis is performed by linearizing about the vehicle equations of motion at 
constant time. This procedure results in a set of third order equations. The autopilot 
gain constants are calculated, and the system transfer functions are established. These 
equations can be reduced to second order by assuming that the winds produce negligible . 
drift in the flight profile. The analysis demonstrates that essentially the same transfer 
functions result when the vehicle motion is approximated by the second order equations. 
Therefore, for simplicity, the second order equations are used to obtain most of the 
results. Later, an integral gain is added to the autopilot in order to guarantee that the 
vehicle will follow the desired pitch program, even if the autopilot and vehicle parameters 
are not accurately known. It is found that integral gain does not substantially reduce the 
effectiveness of the autopilot. 

The three sets of vehicle and autopilot equations (third order, second order, and 
second order with integral gain) are simulated along with triangular wind profiles of 
varying shape and duration. Results are presented for a typical solid propellant launch 
vehicle consisting of a 260 inch solid motor first stage, the SIVB second stage and the 
extended Voyager payload (ref. 3). Maximum deflection angle requirements are pre- 
sented as a function of wind duration. Typical deflection and vehicle attitude profiles are 
also presented with and without limited thrust vector deflection angle capability. The 
equations and results obtained are generalized to apply to other launch vehicles. 


ANALYSIS AND ASSUMPTIONS 

In order to simplify the analysis and allow a closed form solution, the vehicle equa- 
tions of motion are linearized about the nominal trajectory at constant time. A rigid body 
configuration is also assumed. Therefore, vehicle bending, sloshing, and aeroelastic 
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effects are not considered. The angle of attack sensor and the engine gimbal actuator 
are assumed to have no appreciable dynamics. The results presented later show that 
actuator response requirements for the autopilot studied herein are not significantly 
greater than required for a conventional autopilot designed to trim the vehicle. The flight 
path angle can be calculated in the guidance computer from the vehicle position and 
velocity components, which are obtained by integrating the vehicle equations of 
motion. The dynamics of this loop are also neglected in this report. 



Time, sec 


The results are obtained for triangular wind disturbances, with total wind duration 
t (see above sketch) ranging from 3 to 20 seconds. Symmetrical triangular wind shapes 
are investigated extensively, but unsymmetrical shapes are also considered. For the 
unsymmetrical wind shapes, the disturbance buildup time is varied from 3 to 10 seconds. 
The variables used are defined in the symbol list (appendix A). A variety of real wind 
profiles can be approximated by varying the peak wind velocity, duration, and shape of 
the triangular profile. Also, the synthetic wind profiles developed in a work entitled 
’’Terrestrial Environment (Climatic) Criteria Guidlines for Use in Space Vehicles 
Development, 1964 Revision" by Glenn E. Daniels of the NASA Marshall Space Flight 
Center, can be approximated by this method. 

The reduction in thrust vector deflection (TVD) requirements obtained herein applies 
to wind disturbances only. Other TVD requirements, such as TVD required for thrust 
misalinement, center of gravity offset, pitch program, and vehicle dispersions, are not 
affected by the autopilot design chosen. However, these requirements are generally much 
smaller than the TVD required for winds. Therefore, the overall reduction in TVD 
requirements which will be demonstrated is still very significant. 

The autopilot design criteria are as follows: 

(1) The control system undamped natural frequency and damping ratio are specified. 
The damping ratio is set at 0. 707, which results in minimum overshoot and rapid settling. 

(2) The vehicle attitude should follow the nominal pitch program in the absence of 
disturbances. That is, 0/0 =1 at zero frequency. 

r 
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(3) The maximum deflection angle required to maintain stability when disturbances 
(winds) are present is to be minimized. 

The basic vehicle configuration is presented in figure 1. Many of the variables used 
in this report are defined in this figure. As shown in figure 1, the vehicle has aerody- 
namic forces acting on it. These forces may be assumed to be concentrated at a single 
point, the vehicle center of pressure. Generally, the center of pressure is located above 
the center of gravity, which results in an aerodynamically unstable vehicle. The vehicle 
is stabilized by the control system. 



Figure 1. - Basic vehicle configuration. 


The derivation of the vehicle and control equations for the second order, third order, 
and second order plus integral gain systems is presented in appendix B. The equations 
are derived in a cartesian coordinate system whose axes are oriented along and perpen- 
dicular to the vehicle flight path. The origin of the coordinate system is located at the 
center of gravity of the vehicle. Although the equations are derived for the pitch plane, 
they are easily modified for the yaw plane, as shown in the equation derivations. 

The linearized equations of motion in the pitch plane (eq. (Bl)) are: 

A6 = j u c 6 + ij. a a 
a = A0 - Ay + y 


Ay = — [(T - F a ) A0 - (T - F a - mg sin y N ) Ay - T6 + F n a] 
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Av = - g COS y N Ay 


w 

y = smy N 

V N 

where y is the wind angle resulting from the wind disturbance, as illustrated in figure 1. 
The wind angle may be calculated from the wind velocity, nominal flight path angle, and 
vehicle relative velocity, as shown in equation (Blc). Also, since a fixed time analysis 
is used, the wind duration may be related to altitude interval by multiplying by the vehicle 
vertical velocity: 


Altitude interval of wind = rv^ sin y^ 

The nominal flight is assumed to be zero angle of attack with no wind. The control equa- 
tion is: 

6 =K pV K 0 Afl-^Afl-K Ay-K a a 

The desired pitch program is satisfied by proper choice of the gain constants as a func- 
tion of K^. The derivation is presented in appendix B, but the major results are re- 
peated here. The required values of the gain constants are: 


^c 



K 


w n + ^ a 


■e K « 


Kr = 
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K = K - 

y a 
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With these gain constants, the vehicle attitude and control responses become: 
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a e _ M c K p 

6 p (s 2 + 2£o> n s + &> 2 ) 



where 


M- 


1 + ) T + F fi - F a - mg sin y N 


M 


c = 


mv- 


N 


The pole location c is set equal to the zero of AO/6 This reduces the Ad/ 6^ re- 
sponse to second order, which will be discussed later. The pitch program response 
AB/ 6p is specified by adjusting so that the control system follows the pitch program 
in the absence of disturbances. Also, the value of K can be chosen to minimize the 
deflection angle necessary to maintain stability when disturbances are present. 

A second order approximation to the equations of motion is also investigated. This 
simplification is desirable since it permits a better understanding of the system. By 
comparing the second and third order equations, the effect of the additional zero and pole 
can be determined. Also, the effect of an integral gain term added to the system can be 
more easily observed by comparison with a second order system. The second order 
approximation is derived by assuming that the flight path angle does not deviate from the 
nominal (Ay = 0). For this case, the control equation is: 


if 
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Figure 2. - Second-order block diagram. Without integral gain: G(s) = 1. With integral gain: K e = K 0 ; 
G(s) - (1 + Kj/s). p 


5 = Vp - K e Ae - k r a « - K c a 

A block diagram of the second order system is presented in figure 2. Appendix B pre- 
sents the second, as well as the third order equations and the response solutions. Com- 
parison of these response equations demonstrates that the second order equations very 
closely approximate the third order system under investigation. It is shown that the 
large-valued poles and zeroes are identical for the two systems, leaving a small pole and 
zero in the deflection angle response equations. These roots are small and do not con- 
tribute appreciably to the response transients. 

A drawback to the control system design presented thus far is that the vehicle param- 
eters are not precisely known and they vary due to winds. Also, the gain constants can- 
not be set exactly. Due to these factors, the pitch response cannot be followed exactly. 
The second order system was therefore run with an integral gain term added to the 
deflection angle equation. By design, the system with integral gain always follows the 
pitch program at low frequencies regardless of vehicle parameter and gain constant vari- 
ation during flight. The control equation for this case is: 

6 = ( l + ?)vV if,) ' KR4,i - K o“ 

The value of the integral gain constant (Kj) was chosen by specifying the third pole of 
the system (0) and equating coefficients. Since larger 0 values result in a faster pitch 
response and also a larger deflection angle requirement, 0 is chosen as small as pos- 
sible to yield a satisfactory pitch response. The value of 0 was selected to be 


7 













0. 3 radian per second. With this value, it is shown that the deflection requirement does 
not increase appreciably above the requirement obtained without integral gain. 


RESULTS AND DISCUSSION 

The vehicle, trajectory, and aerodynamic parameters used are those of a large 
solid launch vehicle consisting of a 260-inch solid first stage, SIVB second stage, and 
extended Voyager payload. This vehicle is discussed in detail in reference 3. The 
vehicle parameters (table I) were chosen to correspond to maximum dynamic pressure 
flight conditions, since maximum deflection angle requirements occur at about this time. 
However, the reduction in deflection requirements is also applicable for other flight con- 
ditions. The control system undamped natural frequency and damping ratio were specified 
to be 1 radian per second and 0. 707, respectively. These values were taken from refer- 
ence 3. 


TABLE I. - 260-INCH SOLID-VOYAGER VEHICLE PARAMETERS 


Parameter 

Value 

Distance from gimbal point to center of gravity, a, m 

23.6 

Distance from center of gravity to center of pressure, b, m 

17.4 

Axial drag, F^, N 

1. 08xl0 6 

Normal force per angle of attack, F N/rad 

2 ** 

Moment of inertia, I, k g-m 

7.92X10 6 
2. 17x10® 

Mass, m, kg 

1. llxlO 6 

Thrust, T, N 

24. 7x10® 

Velocity, v, m/sec 

583.4 

Nominal flight -path angle, y N , deg 

45.0 


The extended Voyager vehicle was chosen as an example since it has the greatest 
deflection angle requirement of solid fuel vehicles under study (ref. 1). The selection of 
this launch vehicle does not limit this analysis to only one vehicle. Appendix C presents 
a generalization of the vehicle equations developed and demonstrates that the analysis is 
applicable to other solid propellant launch vehicles. 

The linearized equations of motion and control equation were programmed and run on 
a digital computer. Wind duration (t) ranged from 3 to 20 seconds. From reference 2, 
the wind duration which has the largest deflection angle requirement for the launch vehicle 
under study is 4 seconds. Therefore, this t value was investigated extensively. The 
wind duration range was extended to 20 seconds because typical durations are about 10 to 
15 seconds. Also, this expanded range allows for application to other vehicles whose 
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wind duration area of interest may be considerably different. Several values of K fl were 
investigated. The program satisfactorily solved for the pitch and disturbance responses 
for both the second (with and without an integral gain term added) and third order approx- 
imations. 

Figure 3 shows the values of the gain constants required to satisfy the desired trans- 
fer function, as a function of K fl , for the second and third order equations. does 
not exist for the second order case; the other constants are identical for the second and 
third order equations. 

Figures 4(a) and (b) present the maximum deflection ratio, 6-^, (6^ is defined as the 
deflection angle divided by the maximum wind angle) against wind duration for symmet- 
rical triangular wind disturbances with total duration from 3 to 20 seconds. The angle of 



Figure 3. - Variation of feedback gain constants as func- 
tions of angie-of-attack gain constant. 
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attack gain constant (K a ) is varied between 0 and 0. 5. The results in figure 4(a) are for 
the second order equations, while figure 4(b) presents the results for the third order 
equations. The results are similar with the third order system requiring a smaller 
6 r max for small K fl . This is due to drift, which reduces the angle of attack and, 
hence, the deflection angle required to control the vehicle. For both the second and third 
order systems, conventional autopilots (K^ = 0) require significantly larger deflection 
angles than those required for load relief autopilots for wind durations greater than about 
3 seconds. The largest reduction in deflection requirements is obtained by using the 
largest possible value of K^. However, a limiting factor for choosing the magnitude of 
K a is the possible high frequency noise in the angle of attack sensor. Therefore, 

K a = 0. 5 was chosen to be a particular interest since this is the smallest value which is 
large enough to result in a substantial reduction in deflection angle requirements. 

For short wind duration, the = 0. 5 autopilot commands larger deflection angles 
than the conventional autopilot. However, it will be demonstrated in later figures that 
the commanded deflection angle for short wind durations is unnecessarily large (for 
K = 0. 5) and that the vehicle remains stable when a maximum deflection angle (6t, ) 

Cl flj III ct A. 

smaller than the commanded value is used. It will also be shown that the deflection re- 
quirements can be reduced for K a = 0 and large wind durations, but this reduction is 
much smaller than that which is obtained by using a K a = 0. 5 autopilot. 

Figure 5 illustrates the wind angle, angle of attack, and deflection profiles for a 
typical wind with a peak at 3 seconds and a total duration of 12 seconds, for = 0. 5. 



Figure 5. - Wind angle, angle of attack, and deflection angle as functions of time for a triangular wind disturb 
ance. 
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The load relief action of the autopilot is evidenced by the fact that the vehicle angle of 
attack is smaller than the wind angle. For a conventional autopilot, the angle of attack 
and wind angle would be identical for the second order approximation. 

Figures 6(a) and (b) present the 6^ response transients for a symmetrical tri- 
angular wind of 4 seconds duration with K ffl values of 0 and 0. 5. Here, figure 6(a) is 
for the second order analysis and 6(b) for the third order one. It is shown that the higher 
value of commands a larger 6^ for about the first two seconds, which results in a 
lower maximum deflection requirement (6^ = 0.23 for = 0. 5 and 6 R = 0. 31 for 
K = 0). Also, it is interesting to note that the integral of the absolute value of 6^ is 
about 20 percent less for K a = 0. 5 than = 0. This is of particular importance when 
liquid injection is used for obtaining thrust vector control since 20 percent less fuel is 
required for K ff = 0. 5. Also, it can be observed from figure 6 that engine gimbal actua- 
tor requirements have notbeen significantly increased by the use of the load relief autopilot. 

Figure 7 presents 6^ for a symmetrical triangular wind of 4 seconds duration and 
= 0. Two curves are illustrated; one with no limit on Sj^ and the other with a limit 
of 6j^ max = 0. 240. This graph demonstrates that a stable system can be maintained 
when the deflection angle is limited for the = 0 case. The maximum deflection re- 
quirement is 0.31, but 6j^ can be limited to 0.240 and the system remains stable. This 
is accomplished by using 6^ = <5^ max from about 2 seconds to 6. 8 seconds. Prior to 2 
seconds and after 6. 8 seconds, 6^ is less than 6^ max - Trim conditions are essentially 



Figure 7. - Deflection response for a symmetrical triangular wind of 4 seconds, 
with and without deflection -angle limit. Angle of attack gain constant, K Q , 0. 
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Ratio of deflection angle to maximum wind angle, 





(c) Wind duration, 5 seconds. 

Figure 8. - Deflection response with various deflection -angle limits. Angle of attack gain constant, K Q , 0. 5; symmetrical triangular 
wind disturbance. 
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reestablished at about 13 seconds or about 9 seconds after the wind has subsided. There- 
fore, for = 0, the requirement cannot be reduced below about 0. 24. 

Figures 8(a), (b), and (c) present the deflection ratio response for symmetrical tri- 
angular winds of 3, 4, and 5 seconds duration, respectively, with = 0. 5. The deflec- 
tion ratio responses are shown for various values of 6 R max - The 5 R rnay values pre- 
sented on these figures are the smallest values investigated which result in a stable sys- 
tem, plus three additional values. It can be seen from the figure that for t = 3 seconds, 
6 r max can be reduced to 0. 115; for r = 4 seconds, 5 R max = 0. 12; for t = 5 seconds, 
6 tj = 0. 12, and the vehicle remains stable. 

iii ' 

From figures 8(a), (b), and (c), 6 R max = 0. 120 is the minimum value which results in a 
stable system for a wind duration range of 3 to 5 seconds. This value of 5 0 was 

investigated in figure 9 for 6, 8, and 10 seconds duration winds with = 0. 5. Longer 
wind durations were not examined since they never have a commanded deflection require- 
ment greater than 0. 120 (fig. 4(a)). Figure 9 demonstrates that a limit as small as 0. 120 
can be used and still result in stable system. 

Figure 10 presents the ratio of attitude deviation to maximum wind angle (0 R ) re- 
sponse curves for a triangular wind of 4 seconds duration with K ffl values of f) and 0. 5. 
The 0 R response curves for 3, 4, and 5 seconds duration triangular winds are plotted 



Figure 9. - Deflection response for symmetrical triangular winds of 6, 8, and 10 seconds. 
Deflection ratio limit, 0.12; angle of attack gain constant, K Q , 0.5. 
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Angle of attack 



Figure 10. - Attitude response for symmetrical triangular wind disturbance of 
4 seconds. 


in figures 11(a), (b), and (c), respectively. Here, = 0. 5 and curves are presented 
where the maximum deflection ratio is limited to the minimum values found for a stable 
system in the corresponding deflection angle analysis (fig, 8). These curves demonstrate 
that the attitude deviations are acceptable for the minimum 6^ max cases. 

It has been shown that a deflection ratio limit of 0. 120 can be used and still result in 
a stable system for symmetrical triangular winds. This is significant since the require- 
ment for a conventional autopilot designed to trim the vehicle (K^ = 0), is 0. 240. This 
results in a reduction in deflection requirements of about 50 percent for the autopilot 
analyzed. According to reference 2, the deflection requirement can be reduced to 38 per- 
cent of the trim value with an ideal autopilot design. With a more reasonable design, how- 
ever, the requirement was 56 percent of the trim value. Therefore, the closed loop auto- 
pilot analyzed herein has achieved the reduction in deflection requirements indicated by 
reference 2. 

Since actual flight winds are not always symmetric in shape, deflection requirements 
were also studied for unsymmetrical triangular winds with the = 0. 5 autopilot. A 
range of time of peak wind (r 1 ) of 3 to 10 seconds was investigated, with total wind dura- 
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Time, t, sec 

(a) Wind duration, 3 seconds. (b) Wind duration, 4 seconds. 



tc) Wind duration, 5 seconds. 

Figure 11. - Attitude response with and without deflection -angle limit. Angle of attack gain constant, K a , 0.5; symmetrical triangular winds. 


tion ranging from (t^ + 2) to (r^ + 10) seconds. Since the peak deflection requirement for 
= 0. 5 always occurs at the peak of the wind (e.g. , fig. 5), the maximum deflection 
ratio for an unsymmetrical wind with time of peak wind t-^ is the same as for a sym- 
metrical wind of total duration r = 2 tj. Therefore, the maximum deflection ratio 
against wind duration results can be obtained from figures 4(a) and (b). 

The deflection ratio limit of 0. 120 was then tested for each unsymmetrical wind which 
has a commanded deflection ratio greater than this value. It was found that the vehicle 
is stable with 6^ m = 0. 120 for all winds with not less than 4 seconds. For a 
of three seconds, a 6 R of 0. 130 was required for stability. Therefore, the pre- 
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viously obtained reduction in deflection requirements has been slightly degraded for this 
case. However, a of three seconds corresponds to an extremely rapid (and unlikely) 
wind buildup rate. For example, for a peak wind velocity of 60 meters per second, a 
of three seconds corresponds to a wind shear of 0. 049 meters per second per meter. 

That is, the wind velocity increases from 0 to 60 meters per second in an altitude interval 
of 1236 meters. 

Figures 12(a) and (b) present deflection ratio responses for two typical unsymmet- 
trical winds, by using second order equations and = 0. 5. The figures show the un- 
limited 6 r profiles, as well as the stable limited cases. 



Figure 13. - Maximum value of ratio of deflection angle to maximum wind angle as function of symmet- 
rical wind duration. Angle of attack gain constant, K„, 0. 5; second-order system with integral gain. 


Figures 13 to 16 illustrate results for the integral gain system (the second order 
system with the integral gain added). In figure 13, maximum deflection angle ratio 
against wind duration is plotted for symmetrical triangular winds with K a = 0. 5 and 
j3 = 0, 0. 1, 0. 2, and 0. 3 radians per second. The wind duration range is from 3 to 
20 seconds. Observe that the addition of the integral gain results in an increased deflec- 
tion requirement for long wind durations. Also, the higher values of /3 result in larger 
deflection angle requirements for wind durations greater than about 6 seconds. The 
increased 5 R is a result of the integral gain forcing the attitude error to zero. This 
increases the angle of attack and, therefore, results in a larger 6 R . 

Figure 14 compares the 6 R response for the second order analysis with that of the 
third order one which results from the addition of integral gain. It is seen that these 





transients are reasonably similar with a slightly larger deflection requirement for the 
second order system. 

Figure 15 presents the deflection ratio response for the system with integral gain. 
With a wind duration of 4 seconds and = 0. 5, the system was investigated with a 
limited 6^ in order to make a comparison with the corresponding second order analysis. 
The same deflection limits are used as for the second order case; and it is seen that the 
system with integral gain is also stable for these values. 

Since longer wind durations result in an increased deflection angle requirement, a 
wind duration range of 10 to 22 seconds is investigated in figure 16. A stability map for 
the integral gain system is plotted in this figure. Analyzing the results, it is seen that 
for wind durations of up to 20 seconds, a minimum 6 max of 0. 170 is required as com- 
pared to 0. 120 for the second order system. However, for long duration winds, the 
flight parameters are not constant as is assumed by this analysis. Since the analysis is 
made at the maximum dynamic pressure point, 6 required other than at this point 
will be lower for the same wind. Therefore, the 0. 170 value is unrealistically high. In 
addition, /3 can be reduced in the high wind region so as to decrease the angle of attack 
and therefore the deflection requirement. 



Figure 16. - Deflection ratio required for stability. Symmetrical triangular winds; second-order system 
with integral gain; additional pole location, p, 0.3 second"!; angle of attack gain constant, K Q , 0.5. 


CONCLUDING REMARKS 

An autopilot has been analyzed which reduces the deflection angle required for flight 
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through design winds. This reduction is 50 percent of the requirement necessary to fly 
a vehicle under trimmed conditions. The autopilot analyzed requires the feedback of 
angle of attack and flight path angle in addition to attitude error and attitude rate. 

The vehicle parameters simulated are those of a large solid launch vehicle consisting 
of a 260-inch solid motor first stage, SIVB second stage, and an extended Voyager pay- 
load. However, the equations developed may be applied to any launch vehicle. In addition, 
it was demonstrated that the results also apply to other solid propellant launch vehicles. 

Third order equations were used to simulate the system. Also, integral gain was 
investigated. It was shown that the results of both of the above representations were the 
same as approximated by a system of second order equations. 

The 50 percent reduction in deflection requirements obtained applies to the require- 
ment due to wind disturbances only, and not to requirements due to thrust misalinement, 
center of gravity offset, other vehicle dispersions, or pitchover. However, these re- 
quirements are usually much smaller than the wind requirement. For example, for the 
Voyager vehicle, the wind requirement is 2. 3°, the misalinement requirement is 0.25°, 
and the requirement due to the pitch program is 0. 6° (ref. 1). 

A rigid body configuration, as well as perfect angle of attack and flight path angle 
sensors, has been assumed. Further analysis is required to determine the effect of an 
imperfect sensor, actuator dynamics, and vehicle bending and sloshing effects. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 17, 1968, 

125-17-05-01-22. 
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APPENDIX A 


a 


b 

CG 

CP 

c 


GP 

G(s) 

g 


K t 


*P 

k r 


K r 

K e 

K l» K 2 


m 

s 

T 

t 


SYMBOLS 


distance from gimbal point to 
center of gravity, m 

distance from center of gravity 
to center of pressure, m 

center of gravity 

center of pressure 

additional pole in third order 
analysis, sec - * 

axial drag, N 

normal force per angle of 
attack, N/rad 

gimbal point 

gain function 

gravitational acceleration, 
m/ sec^ 

o 

moment of inertia, kg-m 

integral gain constant 

pitch gain constant 

attitude rate gain constant, sec 

angle of attack gain constant 

flight path angle gain constant 

attitude gain constant 

generalized parameters 
(appendix C) 

mass, kg 

Laplace operator, sec - * 
thrust, N 
time, sec 


v velocity, m/sec 

y wind angle, rad 

a vehicle angle of attack, rad 

0 pole which results from the use 

of integral gain, sec - * 

y flight path angle, rad 

A( ) linearized variable 

6 thrust vector deflection angle 

^trim ratio of maximum deflection 

angle to maximum wind angle 
for trimmed flight 

£ damping ratio 

6 vehicle attitude, rad 

0p commanded vehicle pitch at- 

titude, rad 

A dummy variable (appendix C), 

sec - * 

X-^,X 2 generalization variables (appen- 

dix C), sec - * 

H vehicle control parameter, 

c -2 

sec 

ja vehicle aerodynamic parameter, 

01 -2 

sec 

t total wind duration, sec 

wind buildup time, sec 

u> undamped natural frequency, 

rad/ sec 

Subscripts: 

max maximum value 
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N 


nominal 


w 


wind 


R ratio of variable to maximum wind 
angle 

rel relative 


Superscripts: 

time derivative 
second time derivative 
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APPENDIX B 


EQUATIONS OF MOTION 
Third Order Analysis 

The vehicle equations of motion for a rigid body configuration in the pitch plane 
(see fig. 1) are: 


10 = Ta sin 6 + F n ba! 


a = 0 - y + y 


mvy = T sin (0 - y - 6) - F A sin (0 - y) + F n O! cos (0 - y) - mg cos y 
mv = T cos (0 - y - 6) - F A cos (0 - y) - F n a sin (0 - y) - mg sin y 

. v w . 
sin y = — sin y N 

V N 

Although these equations are written for the pitch plane, the analysis can be applied to the 
yaw plane by setting g = 0. 

The nominal trajectory is assumed to be zero angle of attack and zero wind 
(q! n = 0 , y = 0 ). Therefore, 0 N = y N . The nominal deflection angle (5^ is also assumed 
to be zero. Linearizing about the nominal trajectory at constant time, the above equa- 
tions reduce to: 


or 


I A0 = Ta6 + F n ba 


A0 = m c 6 + ii a a 
a = A0 - Ay + y 


(Bla) 

(Bib) 


Ay = — L_ [(T - F a ) A0 - (T - F a - mg sin y N ) Ay - T6 + F n al + JL. cos y N Av (Blc) 

mV N 4 
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Av = - g COS y N Ay 


(Bid) 


y = — sin y N 
V N 


(Ble) 


where 


^ c 



M 


a 


V 

I 


These equations can be applied to the yaw plane by setting g = 0 and y N = 90 degrees. 

From reference 2, Ay is small, which results in a small Av. Also, the variation 
of v leads to a low frequency ' 'phugoid’ * mode which contributes negligibly to the vehicle 
response. Therefore, Av and Av are assumed to be zero in the analysis that follows. 

The autopilot control law utilizes the feedback of all state variables as well as angle 
of attack and flight path angle, as was discussed earlier. The deflection angle is there- 
fore expressed as: 


6 = Vp - K e A0 - k r a0 - K « a - K y A ^ 


(B2) 


Rewriting equations (Bla) and (Blc) by substituting for a, using Laplace notation 
and assuming all initial conditions to be zero: 


(s - ii a ) A 0 + ii a Ay - m c 6 = P a y 
T + F n ' F A " mg sin 


/ T + F n" F A' 


A 6 + I s + 


fflV 


mv. 


N / \ “'N 

( k rS + K e + K ffl ) A 0 + (K - K a ) Ay + 6 = K p 0 p - K^y 



(B2a) 

F 

T 6 = " y 

(B2b) 

mv N mv N 



(B2c) 


Writing equations (B2) in matrix form: 
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s 2 -,, 


T + F n - F A 


K r s + K, ♦ K a 


T + F n - F A - mg sin y N T 


K . K 
y Ot 


V P - K «y 


Taking the determinant of the characteristic matrix, the characteristic polynomial 
(denominator of the transfer functions) is: 

.3 [ T + F „ - f a - m e sil r N , a y T „ 1.2 
L mv N mv N r J 


M T / v T + F - F A - mg sin y N 

+ -^ + — K R + ^ c ( K e + K a ) + ^c — ~ k r s 

mv N ° mv N _ 

+ ^£L(mg sin rN+ (K e + K y )T) + ^-((T + F n -F A )(K r + K 0 ) - mg sin y N (K 0 + K ff )) 
mv N mv N 

Applying Cramer's rule and using the superposition principle to solve for the variables 
A& and 6: 


— numerator: 


— numerator: 


u K_ s + 
^c p 


1 + ] T + F n - F A - mg sin y N 


K - ^c K a> s + — T [ (1 - V T - F A - sin r N ] 
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- — [K 0 fT - F A - mg sin r „) + F n K ] 
mv N *- 


— numerator: 


y N 0 2 0 M a gsiny N 
— s -fi a s + - 

V N 


— numerator: 

y 


k p 


_3 

s + 


T + F n - F a - mg sm 


mv 


N 


- K c s ' 




K qj( t - F A - mg sin y N ) + F n K 


of R 


7 


mv 


N 


-M 


a 


K 


K , + 


R 


mv 


(T - F a - mg sin y N ) 


N 


RT - F A - mg sin y N )K + (T - F A )K ] 

mv N 

Two of the poles of the desired characteristic polynomial are specified by the choice 
of the undamped natural frequency and damping ratio. The third pole is chosen to cancel 
the zero of A0/0 p . This results in nearly identical second and third order transfer func- 
tions, as will be shown later. The desired characteristic polynomial is: 

(s 2 + 2gco n s + co 2 )(s + c) = s 3 + (2£w n + c)s 2 + (2£co n c + co 2 )s + ecu 2 


where 


c = 


1 


mv 


N 



T + F n - F a - mg sin y N 


The feedback gain constants K R , Kg, and are chosen in terms of to satisfy the 
desired characteristic polynomial. Also, the gain constant Kp is chosen so that the 
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control system follows the pitch program in the absence of disturbances. This procedure 
results in: 


K 


R 


2 £w n 


M, 


it w n + ^a v 

K e a 

Mo 


K = K 


a 


Vja 

^c 


K_ 


co 


2 

n 


^c 


With the above values of c and the feedback gain constants, the response equations can 
now be written: 


A£ 


^C K P 


M, 


S + 


1 + T + F - F* - mg sin y N 

V ^c/ 


mv 


N 


(s 2 + 2^w n s + u> 2 )(s + c) 


or 


A 6 _ M c K p 

q 2 2 

p s + 2tw n s + w n 


M 


1 + j T + F n - F A - mg sin y N | 


(Ma M c K a ) 


s + 


M, 


A£ 

y 


mv 


N 


(s 2 + 2£u> n s + co 2 )(s + c) 
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or 


A 6 _ ^a ^c^a 
y s 2 + 2£a> n s + CO 2 


K , s 3 + T + F n- F A- m S sin ^N„2 
6 P l mV N 


s -l‘a s + ^a 


g Sin y N > 


N 


e. 


(s 2 + 2£a> n s + w 2 )(s + c) 


or 


^ /J \ /I * sin ?n\ 

(s 2 + 2^a> n s + a> 2 )(s + c) 


y (s 2 + 2fco n s + a> 2 )(s + c) 


. < -K^s' 




o 

2 ?w n + 

Mn 


V T - f A - mg sin r N ) - 


mv 


N 


- H 


a 


w n + ^ 


i 2£a; ri 

- K a - — l — (T - F a mg sin y N ) “ 


M, 


mv 


N 




M 


a 


mv 


r A .2 

(T - F a - mg sin y N ) 


N 


VA£ - K a ) + (T - f a )(k ct - 12 .) 




M, 


•A 
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y (s 2 + 2£u> n s + 0 > 2 )(s + c) 


M* + 2Jw„s + ^ (<o 2 + M a - 

^ r r 1 r 


mv 


< s + 


N 


(-mg sin y N ) + -2. (T - F A - mg sin y^ 


w n + 
M« 


- K 


a 


Second Order Analysis 

The linearized equations of motion (eqs. (Bla), (Bib), and (Blc)) presented in the 
first section of this appendix may be approximated by assuming that the flight profile does 
not deviate from the nominal profile, due to the presence of winds; that is, Ay = 0. This 
is a valid assumption for most large launch vehicles, as is shown in reference 2. The 
resulting linearized equations of motion are second order and are given by: 

A6> = n c 6 + n a a 


a = A0 + y 


Vp - K e Ae 


K R A6 - K a a 


Writing these equations in matrix form and solving for the A 9 and 6 responses as in 
the third order case results in: 


A 0 ^*5? “n 

% s 2 + ii c K R a + n c (K e + K a )- n a s 2 + 2fo) n s + co 2 

A0_ M ffl -M c K a 

v 2 2 

y s + 2£a> n s + o> n 
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I 






0 p s 2 + 2fco„s + co 2 
b n n 


6 

y 


K. 


s 2 + — 2£a> s+ — (a) 2 + jut)-jLt K 
a; T ^ n V, n ^ay ^ar a 


M a ( . .2 


s 2 + 2fa> n s + a> 2 


where and Kg are chosen to satisfy the desired natural frequency and damping ratio 
and Kp is chosen to satisfy the desired pitch response. 


K R = 


2 £co 


n . 


M, 


K 0 = ' 


w n + ^a 


Ct> 


- K ; K = — H 
a ’ n 




P 




Comparison of Formulations 

Comparing the results of the second and third order analyses, it is seen that the 
feedback gain constants, Kp, K^, and Kg are identical for the two formulations. An 
additional gain constant, K^, exists for the third order formulation (and also in the real 
case). By setting the additional pole of the third order system equal to the zero of the 
pitch response, A 0/0 and A0/y are also identical for the two systems. The small 

r 

zero and small pole occurring in the third order system do not contribute measureably 
to the deflection angle response. The vehicle response to wind disturbances can therefore 
be accurately simulated by using second order equations. 


Integral Gain Analysis 

In order to more closely approximate an actual control system, an integral term is 
added to the control variable of the second order system. This use of an integral gain 
term in the autopilot guarantees that the vehicle will follow the desired pitch program 
even if the autopilot or vehicle parameters are not accurately known. With 

6 = Vp - K e Ae - K R s Ae - K a a 
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let 


V k s = ( 1 + tH 


Then the deflection angle becomes 


6 = ^ Kp(e p - A 9) - K r s A9 - K a (A8 + y) 


s6 = (s + Kj)Kp(e p - A 6) - K r s 2 Ad - K a (A 9 + y)s 
where a = A6 + y was used. Also, 


A 9 = |U C 6 + ii a a 


or 


(s 2 - n a )A 9 - |i c 6 = M a y 

Rewriting the above equations and applying Cramer's rule to solve for the gain constants 
and transfer functions results in: 

|(s + Kj)K p + K r s 2 + K^s] A0 + s6 = (s + Kj)K p (? p - K^sy 
(s 2 - n a )A 8 - M c 5 = M a y 

A9_ _ __ M c (s + K l)Kp 

0 p s 3 + !i c K R s 2 + [^ c (K p + K a ) - fxjs + fXcKjKp 

The denominator must be set equal to 

(s 2 + 2£u> n s + cu^(s + j8) = s 3 + (2£u> n + /3)s 2 + (2fo> n /3 + w 2 )s + 0o) 2 

Therefore, 
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2tu 8 + o> + a - a K 
fc rr n ' a ‘ c a 


K R “ 


2£w n + 8 


A u K 

p 

Discretion must be used in the selection of /3. Larger 8 values result in both a faster 
pitch response and a higher deflection angle requirement for flying through winds. There- 
fore, 8 is assigned the smallest value possible which yields a satisfactory pitch response. 
In this case, 8 = 0. 3 radian per second was chosen. The deflection requirement does 
not increase appreciably with this value. 

The vehicle attitude response to a wind disturbance becomes: 

AS K ~ K a )s 


y (s 2 + 2£co n s + cu 2 )(s + 8) 


and the deflection angle transfer functions are: 


2 2 

3 2^rv/ 2 \ ^ rv * 

H 2 « S ' 3 + s 2 - feu fi + c 2 + n - M K )s - ,3a, 

M c ^c ^c ^c 

(s 2 + 2^co n s + a, 2 )(s + 8) 

- K s 3 tt°L (2 + /3)s 2 + --9L feu 8 + + M C K )s + ^ ,3a, 2 

5 _ ^c ^c 11 c _ 

y ( s 2 + 2fco n s + o, 2 )(s + 0) 
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APPENDIX C 


GENERALIZATION OF AUTOPILOT ANALYSIS 

By applying the second order approximation of appendix B, the deflection angle to 
wind response equation is a function of only two vehicle parameters, /i^ and p. c , and the 
autopilot natural frequency, co n . 



M 


co„ 




Kq! “ + • s + \i a I — + 






M, 


2 2 
s + 2£co n s + w n 



where it is assumed that the same damping ratio is desired for all vehicles. 

Two constants, and Kg, will now be defined in terms of the three parameters 
mentioned above. Letting 


co 


2 

n 


= K 1^ 


K 


Of 



^_a 

^c 


and Kg can be evaluated. For the 260-inch solid Voyager vehicle with K^ = 0. 5, 
Kj = 1. 57 and Kg = 2. 1. Substituting Kj and Kg into the (5/y)(s) equation results in 
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Defining (S trim /y) = -(m 0 /m c ) 



For a second vehicle, define y/i a - Xg and assume a wind profile with the same shape 
but scaled in time so that the wind duration is Then 



From the definition of a Laplace transform, if 

y{f(t)} = F(s) 

then 
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J 5 f{f(xt)} = AfQ 


Therefore, 



^ 1 ^ 1 ^”) 


Xff (\2^) 


That is, when using a different vehicle (vehicle 2), the wind duration (r) need only be 
transformed by multiplying by (A-j/Ag) and the appropriate deflection ratio can be read 
from the figures presented. 

In order for this generalization to be valid, = 1. 57 must result in an acceptable 
u> n (so not to excite the vehicle bending modes). Table n illustrates the |i a and p, c 
values for two other typical solid propellant rocket vehicles, the 260-inch solid. Apollo 
and SSOPM vehicles (ref. 1), in addition to the 260-inch solid-Voyager vehicle. The 
260-inch solid-Apollo and Voyager vehicles are identical except for different payloads. 
The SSOPM vehicle has solid first and second stages and is designed to deliver one 
million pounds (450 000 kg) to a low earth orbit. The natural frequency for each vehicle 
is calculated and is seen to be reasonable (not appreciably different from the value 
selected for the Voyager vehicle). Any other vehicle can also be tested by solv- 
ing for the u> n value by using the generalization parameters, and Kg, and the 


TABLE n. - VEHICLE PARAMETERS FOR GENERALIZATION 
OF AUTOPILOT DESIGN 


Parameter 


Vehicle 



Voyager 

Apollo 

SSOPM 

Vehicle aerodynamic parameter, 
M a , sec -2 

0.636 

0.436 

0. 171 

Vehicle control parameter, p. , 

-2 c i 

sec 

2.67 

2.67 

1.4 

Ratio of aerodynamic parameter to ' 
control parameter, 

0. 238 

0. 163 

0. 122 

Undamped natural frequency, o> n , 
rad/ sec j 

1.0 

0. 83 

0. 52 
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vehicle parameters and p c . K the calculated a> n is acceptable, the user merely 
needs to translate the abscissa of the figures. If the u> n value is not acceptable, new 
and Kg values can be solved for and new deflection angle to wind duration curves 
must be plotted. Based on the results presented in reference 2, it is expected that the 
same 50 percent saving in deflection requirements would still be obtained. 
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